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Introduction 


The  goal  of  this  project  is  to  perfor  m  initial  development  of  a  new  sensor  syste  m  for 
ensuring  proper  long  term  fitting  of  a  lower  limb  prosthesis  by  providing  real-time  monitoring  of 
pressure  distribution  at  the  body-prosthesis  interface.  To  date,  we  have  completed  t  he  design, 
fabrication  and  optimization  of  the  base  sensi  ng  platform,  the  sensor  ,  and  currently  have  a 
paper  submitted  for  review  and  possible  publication.  Specifically,  we  have  constructed  a 
sensing  layer  capable  of  monitoring  force  on  hard  surfaces  (see  Appendix  A  for  the  submitted 
article).  To  evaluate  the  efficacy  of  the  sensor,  we  constructed  an  automated  pneumatic 
mechanical  loader,  implemented  the  interface  controlling  circuits  for  the  loader  and  developed 
an  algorithm  for  data  analysis.  Additionally,  a  detection  circuit  has  also  been  designed  and  i  s 
being  finalized.  The  current  phase  of  the  project  involves  advancing  the  base  technology  into  a 
grid-like  mapping  se  nsor,  constructing  the  excitation  and  detection  circuits,  modifying  and 
optimizing  the  algorithm  and  embedding  the  sensor. 
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Body 


The  focus  of  this  project  is  the  d  evelopment  of  a  prototype  sensor  system  capable  of 
measuring  pressure  pr  ofiles  at  the  socket-stump  interface  of  a  lower  limb  prosthesis.  T  he 
specific  aims  of  this  pro  ject  are  (1)  Fabrication  of  the  Sensing  Layer,  (2)  Develop  ment  of  the 
Detection  System  and  (3)  System  Performance  Evaluation. 

Specific  Aim  1 :  Fabrication  of  the  Sensing  Laver 
Task  1A:  Sensing  Laver  Fabrication 

We  previously  reported  on  three  investigated  techniques  for  fabricating  the  sensing  layer 
and  at  this  phase  of  the  project  the  use  of  Me  tglas  is  still  ideal  due  to  its  high  stre  ss  sensitivity; 
however,  a  highly  desirable  technique  for  fashioning  sensors  from  Metglas,  screen  printing  and 
etching,  has  been  tested.  Additionally,  this  technique  has  also  allowed  for  preliminary  testing  of 
non-uniform  shaped  sensing  strips  capable  of  measuring  four  loading  points.  Furthermore,  initial 
investigation  into  embedding  the  Metglas  sensors  has  also  been  performed. 

Screen  Printing  and  Etching:  Previously,  a  sh  earing  method  was  utilized  to  fashion  sensin  g 
strips.  While  the  shearing  technique  produces  se  nsors  with  good  magnetic  and  magn  etoelastic 
properties,  it  is  difficult  to  control  the  reproducibility  thus  creating  a  challenge  for  making  a 
uniform  sensing  layer.  As  a  result,  a  new  technique  for  fabricating  these  strips  based  on  screen 
printing  and  etching  has  been  developed. 

The  screen  printing  and  etching  technique  first  involves  the  fabrication  of  an  ena  mel 
mask,  with  the  desired  designs,  that  is  placed  over  a  fine  wire  mesh  screen.  Once  the  screen 
and  mask  a  re  prepared,  the  strips  are  adhered  to  a  p  olycarbonate  block  using  double  sided 
adhesive  film,  and  the  n  placed  on  the  chuck  of  the  screen  printing  machine.  Tests  are  then 
performed  to  ensure  t  hat  the  ap  propriate  amount  of  pressure  is  applied  by  each  of  the 
squeegees,  with  the  first  squeegee  being  used  for  spreading  applied  paint  and  the  second  being 
used  for  pushing  applied  paint  through  the  screen  onto  the  strips.  Following  this,  a  mixture  of 
5900  Series  Enamel  Plus  Gloss  Screen  In  k  (Nazdar  Shawnee  Facility,  8501  Hedge  La  ne 
Terrace  Shawnee,  KS  U.S.A.)  and  9050  Retarder  Thinner  (Nazdar  Shawnee  Facility,  8501 
Hedge  Lane  Terrace  Shawnee,  KS  U.S.A.)  is  applied  on  the  screen,  around  the  designs  and  on 
the  squeegees.  The  machine  then  takes  over  with  an  automated  paint  application  pr  ocess.  The 
strips  are  then  baked  at  85  °C  fo  r  one  hour  and  allowed  to  coo  I.  The  strip  now  has  the 
appropriate  protective  coatings  and  can  be  etched. 

Prior  to  etching,  a  small  beaker  is  filled  with  ferric  chloride  and  placed  onto  a  hotplate  set 
to  45  °C.  While  the  solution  is  h  eating  up,  the  screen  printed  strips  are  sheared  in  order  to 
produce  smaller  str  ips  each  containing  one  coated  design.  These  smaller  strips  are  then 
attached  to  a  specially  designed  a  pparatus  for  holding  the  strips  durin  g  the  etching  process 
using  double  sided  adhesive  tape.  The  apparatus  consists  of  a  polycarbonate  block  with  a  hole 
drilled  through  the  bottom  in  order  to  accommo  date  a  magnetic  stir  stick.  This  is  n  ecessary  in 
order  to  not  only  provi  de  constant  motion  of  the  strips  in  the  etching  solution,  but  to  also 
expedite  the  process.  Once  the  solution  reach  es  45  °C,  t  he  apparatus  is  submerged  and  th  e 
magnetic  stirring  is  turned  on.  It  is  important  to  observe  the  etching  process  every  10  minutes, 
especially  near  the  20  minute  mark,  because,  while  it  is  necessary  to  ensure  that  all  uncoate  d 
material  is  removed,  if  the  strips  are  etched  for  too  long  the  exposed  sides  of  the  coated  areas 
can  be  over  etched  resulting  in  curving  and  other  undesirable  shape  effects  at  the  edges. 

Non-Uniform  Sensor  Desiqn  :  The  use  of  screen  printin  g  and  etching  has  a  Iso  inspired  the 
fabrication  and  preliminary  testing  of  sensing  strips  with  non-uniform  geo  metry  between  th  e 
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front  and  back  halves  of  the  strip  (see  Fig.  1).  In  theory,  when  the  length  of  the  sensing  strip 
moves  away  from  a  detection  coil,  the  strength  of  the  captured  magnetic  field  from  a  particular 
portion  of  the  strip  decreases. 


Fig.1.  Non-symmetric  sensing  strip  design  capable  of  acting  as  a  2x1  sensor,  with  regions  1  and  2 

denoting  different  force  sensing  areas. 

As  a  result,  the  back  hal  f  of  a  strip  will  appear  to  have  a  weaker  overall  response  than 
the  front  half.  This  mea  ns  that  the  front  and  back  halves  of  a  strip  can  be  thought  of  as  two 
different  sensing  elements,  theoretically  allowing  a  sin  gle  strip  to  measure  f  orces  at  two 
pressure  points. 

In  order  for  the  sensing  strip  to  function,  there  must  be  a  distinguishable  difference  in 
sensor  response  when  applying  loads  from  the  front/back  h  alves  of  the  strip  (i.e.  region  1  and 
region  2  respectively).  While  this  effect  has  been  demonstrated  in  rectangular  strips,  the  result 
of  loading  only  region  2  is  not  sensitive  enough  for  this  a  pplication;  however,  if  the  strip  is 
rotated  180  degrees,  then  region  1  becomes  region  2  and,  by  quantifying  the  effect  of  loading  in 
these  orientations,  force  can  be  monitored  in  both  regions  (in  actual  use  there  will  be  detection 
coils  monitoring  both  ends  of  the  strips).  Unfortunately,  the  geometric  symmetry  of  a  rectangular 
strip  presents  a  problem  in  that  th  e  response  of  the  rotated  strip  is  the  same  as  t  he  unrotated 
response.  Additionally,  even  if  the  loading  can  be  determined,  which  region  the  loading  is  from 
cannot  be  discovered  without  visual  inspection  of  the  test  conditions,  whereas  in  the  ideal  case 
only  analysis  software  will  be  utilized  to  determine  the  force  location. 

To  solve  this  problem,  half  rect  angular-half  trapezoid  shaped  sensing  strip  s  were 
proposed  and  preliminary  testing  of  sheare  d  sensing  strips  has  been  performed.  Fig.  2 
illustrates  the  response  of  a  rectangle  and  a  trapezoidal  strip  when  loaded  in  region  1  and  2,  for 
both  rotated  and  unrotated  orientations,  while  Fig.  3  illustrates  a  rectangular/trapezoidal  sensing 
strip  measured  while  loading  regions  1  and/or  2  under  both  orientations. 


Total  Applied  Load  (lbs) 


Total  Applied  Load  (lbs) 


Fig.  2.  (left)  Preliminary  testing  demonstrating  the  difference  in  sensor  response  when  loading  regions  1 
and  2  of  a  rectangular  and  trapezoidal  sensing  strip  while  measuring  from  the  front  (unrotated)  and  back 
(rotated)  orientations. 

Fig.  3.  (right)  Preliminary  testing  of  a  trapezoidal  sensor  while  loading  regions  1  and/or  2  from  the  relative 
front/back  of  the  strip. 
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As  can  be  seen,  measurements  fro  m  the  rectangular  strip  have  little  difference  in  either 
orientation  when  loading  in  region  1  and/or  2.  On  the  other  hand,  the  trapezoidal  sensor  design 
illustrates  a  distinct  d  ifference  not  only  between  loading  in  regions  1  and  2,  but  also  between 
both  sensor  orientations.  It  is  theorized  that  this  results  not  only  from  a  difference  in  the  amount 
of  material  under  load,  which  contributes  to  the  overall  magnetic  fie  Id  strength,  but  also  th  e 
effect  that  the  non-symmetric  geometry  has  on  t  he  magnetic  field  lines.  This  is  currently  only  a 
theory  and  has  yet  to  be  investigated  further. 

In  conclusion,  the  screen  printing  and  etching  process  allows  for  precise  repeat  able 
production  of  sensing  elements  for  grid  fabrication.  Moreover,  the  process  allows  for  the 
exploration  of  the  effects  of  non-uniform  stri  p  design  on  overall  sensor  performance.  The  next 
step  in  this  process  is  to  fully  characterize  a  two  strip  trapezoidal  sensor  design  for  measuring 
four  pressure  positions. 


Embedded  Sensor:  Preliminary  investigation  s  into  embedd  ing  the  sensors  have  also  been 
performed.  Sensing  strips  were  embedded,  in  soft  polyurethane  rubber,  inside  of  a  well  milled 
into  a  piece  of  polycarbonate.  Fig.  4  displays  results  from  preliminary  embedded  sensor  testing 
and,  unlike  data  collected  when  the  sensor  was  placed  on  a  hard  surface  (see  Fig.  5a-c  of  the 
submitted  journal  article  in  Appendix  A),  there  is  an  initial  downward  trend  in  the  data  followed 
by  an  upswing. 
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Fig.  4.  The  response  of  an  embedded  rectangular  sensing  strip  under  loading 


Although  this  effect  has  not  been  thoroughly  investigated  yet,  it  is  theorized  that  this  may 
be  a  result  of  sensor  movement  during  compression  of  the  embedding  material  and/or  th  e 
effects  of  tensile/compressive  forces  on  the  sensor  during  material  compression.  These  results 
will  be  further  explored  following  the  development  of  the  sensing  grid. 

Task  1B:  Sensing  Strip  Optimization 


The  new  fabrication  technique  with  screen  printing  has  significantly  reduced  the  internal  stress 
of  the  sensing  strip.  As  a  result,  th  e  as-fabricated  strips  h  ave  relatively  good  ma  gnetic  and 
magnetoelastic  properties.  Nevertheless,  to  ensure  all  residue  stress  in  the  strips  was  removed, 
all  sensors  were  annealed  in  vacuum  (-600psi)  at  250  °C  for  1  hr. 

Specific  Aim  2:  Development  of  the  Detection  System 


Task  2A:  Detection/Excitation  Circuit  Fabrication 
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Excitation  Circuit:  To  eliminate  the  need  for  external  AC  and  DC  power  supplies,  a  driving  circuit 
was  designed.  AWein  Bridge  oscillator  gener  ated  a  sinusoid  which  was  then  fed  through  a 
fourth  order  multiple-feedback  filter.  Next,  an  amplifier  with  a  class  AB  output  stage  employing 
TIP42  and  TIP43  high-power  transistors  provided  the  power  amplification  necessary  to  obtain 
field  strengths  comparable  to  the  original  AC  supply  and  audio  amplifier  combination. 

For  the  DC  supply  circuitry,  a  DAC  output  from  the  microcontroller  was  amplified  with  a 
TIP33C  transistor  p  laced  in  the  fe  edback  loop  of  an  op  amp.  Use  of  the  DAC  voltage  was 
desired  to  enable  a  programmable  DC  biasing  sweep  of  the  sensors.  The  LF356  was  selected 
for  coupling  with  the  TIP33C  due  to  its  similari  ty  to  the  LF353  with  increased  capabilities  when 
driving  reactive  loads. 

Harmonic  Detection  Circuit:  A  more  efficient  method  of  detection  wa  s  required  to  combat  the 
problem  of  slow  detection  speeds  incurred  by  sampling  coils  in  series.  To  that  end,  a  method  of 
sampling  in  parallel  was  devised  which  involved  multiple  lock-in  detection  circuits  measured 
concurrently  by  a  microcontroller.  A  dditionally,  since  only  the  second  harmonic  is  desired,  t  he 
detection  frequency  was  set  at  400  Hz  to  eliminate  the  need  for  an  entire  spectrum. 

Lock-in  amplification  involves  the  multiplication  of  the  in  put  signal  with  a  reference 
sinusoid  followed  by  low-pass  filtering.  Multiplication  of  sine  waves  has  the  effect  of  generating 
two  new  sine  waves:  one  at  a  frequency  -  f2  and  one  at  +  f2.  When  the  reference  and  input 
sinusoid  are  of  the  sa  me  frequency,  the  result  is  a  DC  voltage  plus  a  sinusoid  at  the  square  of 
the  signals’  frequencies.  The  result  of  multiplying  all  signals  of  dissimilar  frequency  will  be  an 
AC  signal  component.  By  low-pass  filtering,  the  DC  voltag  e  can  be  isolated,  correlating  to  the 
RMS  voltage  of  the  input  component  at  the  frequency  oft  he  reference.  Because  the  multiplied 
signal  results  are  dependent  on  th  e  phase  be  tween  the  reference  an  d  input  sin  usoids,  an 
additional  step  must  be  taken  to  repeat  the  multiplication  and  filtering  process  with  a  reference 
offset  by  90  degrees  from  the  original.  It  is  then  possible  to  extract  the  magnitude  information  by 
applying  the  Pythagorean  Theorem  to  the  unshifted  and  shifted  results.  Similarly,  the  phase  can 
be  determined  by  taking  the  inverse  tangent  of  the  signal  in  relation  to  the  reference. 

The  input  signal  from  t  he  coil  was  first  acquir  ed  with  an  AD8220  instrument  a  mplifier 
input  stage  for  high  CMMR.  The  instrument  amplifier  was  followed  by  a  60  Hz  notch  filter  and  a 
fourth-order  band-pass  filter  with  a  gain  of  100  at  a  center  frequency  of  400  Hz.  Th  e  amplified 
and  filtered  signal  was  then  multiplied  with  a  reference  signal  generated  by  a  Wein  Bridg  e 
oscillator  circuit  at  400  Hz.  The  multiplier  selected  for  this  application  was  the  AD835  4-quadrant 
multiplier  due  to  its  fast  settling  time  ,  low  noise,  high  input  impedance  and  minimal  requirement 
for  external  components.  The  inpu  t  signal  was  also  multiplied  by  a  version  of  the  reference 
frequency  fed  through  a  90  degree  phase-shifting  circuit.  After  multiplication,  second-order  low- 
pass  Butterworth  filters  extracted  the  DC  signal  from  the  two  outputs.  T  he  current  status  of  the 
circuit  proves  its  ability  to  extract  the  DC  level  based  on  frequency  and  phase  but  has  yet  to  be 
integrated  with  a  microcontroller  and  the  rest  of  the  system. 

Task  2B:  Detection  Coil  Fabrication 


New  detection  coils  and  a  new  faceplate  were  constructed.  A  custom  built  manual  coil 
winder,  with  a  rotary  co  unter,  provided  a  simple  method  for  ensuring  a  consistent  number  of 
turns  on  ea  ch  coil  and,  as  a  result,  the  same  or  nearly  t  he  same  impedance.  Additionally, 
previous  errors  in  detection  coil  construction  were  minimized  by  using  a  VFIash  3D  printer  to 
construct  the  coil  frames.  Wound  coils  were  at  tached  to  a  new  faceplate,  fabricat  ed  using  a 
CNC  Micro  Milling  Machine,  which  was  then  screwed  onto  the  excitation  coil. 

Task  2C:  Alqorithm  Development 
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We  previously  reported  the  development  of  a  linear  multivariable  analysis  techniq  ue  to 
determine  the  applied  load  as  a  fu  notion  of  2"'^  order  harmonic  amplitude.  However,  the  new 
automated  loader  has  allowed  for  the  collection  of  more  precise  data  which  has  more  clearly 
shown  what  the  relationship  between  applied  load  and  2'^^  order  harmonic  amplitude  appears  to 
be. 

Linear  Multivariable  An  alysis:  The  previously  p  resented  algorithm  relied  on  the  u  se  of  linear 
multivariable  analysis.  In  essence,  analysis  of  data  collected  prior  to  the  automated  systems  led 
to  the  conclusion  that  t  he  relationship  between  applied  loa  d  and  2'^'^  order  harmonic  amplitude 
could  be  linearized  and,  thus,  placed  into  a  series  of  matrices  which  could  be  easily  manipulated 
to  find  force.  The  previously  used  model  was  abandoned  when  it  was  discovered  that  the 
relationship  between  applied  force  and  2"'^  order  harmonic  amplitude  more  closely  resembled  a 
decaying  exponential  function  which,  while  capable  of  being  linearized,  has  multiple  roots,  thus 
making  the  linearizatio  n  unreliable.  Additionally,  an  e  xponential  curve  model  fit  s  with  the 
theoretical  description  of  the  permeability  change  in  response  to  force  described  in  section  2.3 
of  the  appended  journal  article  in  Appendix  A. 

Superposition:  The  new  algorithm  relies  upon  the  principle  of  superposition.  Due  to  the  close 
proximity  among  the  magnetoelastic  str  ips,  the  stress  response  of  each  magnetoelastic  strip 
was  found  to  be  dependent  on  the  loading  conditions  of  neighboring  strips,  in  addition  to  its  own 
loading  condition.  Therefore,  assuming  the  strip-to-strip  cross  interference  is  cumulative  a  series 
of  equations  were  devel  oped  to  relate  applied  force  on  a  strip  to  measured  2  order  harmonic 
amplitude  at  each  coil.  A  detailed  description  of  these  equations  and  th  e  theory  behind  them 
can  be  found  in  section  2.4  of  the  appended  journal  article  currently  in  review  (see  Appendix 
A). 

Specific  Aim  3:  System  Performance  Evaluation 

According  to  our  initial  plan.  Specific  Aim  3  will  be  conduct  ed  in  full  after  completion  of 
Specific  Aim  1  and  2.  Specifically,  we  will  e  mploy  the  sensor  fabr  icated  in  Aim  1  into  th  e 
detection  system  developed  in  Aim  2.  The  performance  of  the  whole  technology  will  then  be 
evaluated  and  optimized. 

Even  though  Aim  3  is  not  yet  addressed  completely,  certain  portions  of  the  system  have 
been  reevaluated,  out  of  necessity  ,  and  found  to  be  insufficient  for  t  he  continuation  of  the 
project.  Specifically,  previously,  the  use  of  a  four  shaft  manual  mechanical  loader  was  reported. 
This  loader  proved  to  be  insufficient  for  continuing  work  and  a  new  mechanical  loader  and 
detection  system  have  been  developed.  The  new  system  consists  of  a  fully  aut  omated  four 
piston  pneumatic  mechanical  loader  and  an  a  dditional  automated  three  coil  detect  ion  system, 
with  the  appropriate  circuit  elements  and  computer  control  programs. 

Automated  Mechanical  Loader:  Fig.  5  illustrates  the  new  mechanical  loader.  The  outer  frame  of 
the  loader  was  comprised  of  four  stainless  steel  rods  screwed  into  a  delrin  base  and  an 
aluminum  top.  The  stainless  steel  rods  provided  necessary  structural  support  and  the  aluminum 
top  allowed  for  metal  threading  for  screwing  in  the  pneumatic  piston  s,  which  was  necessary 
since  screwing  metal  into  plastic  or  plastic  into  metal  is  never  ideal  especially  when  the 
plastic/metal  will  be  under  loading.  The  rest  of  the  system  consisted  of  delrin  in  order  to  prevent 
the  construction  of  a  metallic  frame  which  could  create  a  shielding  effect. 

Force  is  applied  by  four  delrin  rods,  with  associated  load  cells,  attached  to  stainless  steel 
cylinders,  screwed  onto  the  piston  shafts  by  slotted  connector  pieces  custom  built  using  aV- 
Flash  FTI  230  Desktop  Printer  (3D  Systems  Corporation  333  3D  Systems  Circle,  Rock  Hill  SC 
29730).  Pressure  build  up  in  the  pistons  causes  the  shafts  and  attached  cylinders  to  extend.  As 
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the  cylinders  move  down,  the  hanging  rods  and  load  cells  also  move  until  contacting  the  sensor. 
As  a  result  of  the  slot  in  the  connector  pieces,  th  e  cylinders  continue  to  move  until  reaching  the 
load  cells,  thus  applying  force.  This  not  only  c  reates  a  true  zero  loading  condition  ,  since  the 
delrin  rods  hang  from  t  he  cylinders  and  do  not  initially  contact  the  sensor,  but  also  allows  the 
application  of  different  force  profiles  in  four  regions.  In  addition,  a  stepper  motor,  not  pictured  in 
Fig.  5,  is  attached  to  the  plate  on  which  the  sensor  sits,  allowing  for  computer  controlled  rotation 
of  sensors  when  testing  from  the  front/back  with  a  single  set  of  detection  coils. 


Fig.  5.  The  mechanical  loader  houses  the  sensor  and  is  capable  of  providing  different  loading  profiles  in 

four  regions. 

Air  flow  to  the  pistons  is  controlled  by  a  series  of  Clippard  EVP  series  proportional 
control  valves  (Clippard  Instrument  Laboratory,  Inc.  7390  Colerain  Ave.,  Cine  innati,  Ohio 
45239).  These  solenoid  type  valves  control  airflow  as  a  function  of  applied  voltage,  with  zero 
volts  completely  closing  the  valve  and  ten  volts  completely  opening  the  valves.  Air  coming  from 
the  main  source  enters  into  the  control  box  where  it  is  split  between  two  T  junctions,  essentially 
creating  four  lines  of  airflow.  Each  of  the  four  air  pathways  i  s  then  split  again  upon  meeting  the 
first  of  two  EVP  valves.  The  first  valve  acts  as  a  main  control  for  overall  air  volume  control  and 
splits  the  air  between  a  piston  and  another  EVP  valve.  The  second  valve  acts  as  the  primary 
control  unit  for  force  application.  By  opening  or  closing  the  second  EVP  valve,  the  total  volume 
of  air  to  the  piston  cha  nges.  A  higher  volume  of  air  flow  to  the  piston  will  result  i  n  greater  air 
pressure  and,  thus,  a  higher  load,  and  vice  versa.  In  total,  the  airflow  control  system  consisted 
of  eight  EVP  valves  connected  to  four  pneumatic  pistons. 

Loader  Control  Circuitry :  The  hard  ware  design  for  then  ew  mechanical  loader  included  four 
parallel  digital  proportional  plus  integral  (PI)  control  units  realized  through  the  integration  of  a 
microcontroller  (MSP430fg4618),  an  array  of  transistors  for  EVP  actuation  (2SD882)  and  load 
cells  (FC231 1 )  for  feedback  on  the  current  loading  state.  Additionally,  offset,  gain  and  filtering 
stages  were  also  implemented  for  design  optimization  and  calibration. 

For  acquisition  of  load  data,  four  Measurement  Specialtie  s  FC23  250  Ibf  Compression 
Load  Cells  placed  in  series  with  the  applied  load  and  actuating  pneumat  ic  cylinders.  Differential 
amplifier  configurations  (LF353)  then  converted  the  amplified  millivolt-level  signals  to  achieve  an 
output  range  of  0  to  2.  5  V,  which  could  be  measured  by  the  microcontroller’s  analog-to-digital 
converter.  Due  to  repeatability  issue  s  with  the  load  cell  offset  voltages,  a  small  offset  (approx. 
40  mV)  wa  s  added  to  the  amplified  load  cell  voltage  to  account  for  its  maximum  negative 
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voltage  offset  and  thus  assure  that  all  microcontroller  inputs  were  positive.  In  combination  with 
the  negative  offset-nulling  circuit,  t  he  microcontroller  software  automatically  compensated  fo  r 
the  new  positive  offset  at  power-up  and,  as  de  sired,  between  tests.  This  feature  significantly 
improved  the  reliability  of  the  force  data,  especially  when  collected  over  long  durations. 

In  order  to  impleme  nt  the  desired  loadin  g  control  algorithm,  a  MSP430f  g4618 
microcontroller  unit  (MCU)  was  utilized.  This  microcontrolle  r  was  selected  due  to  it  s  250  ksps 
12-bit  analog-to-digital  converter  (ADC),  low  power  operation,  16  MHz  clock  fr  equency  and 
integrated  SPI  and  UART  peripherals.  Additionally,  previous  knowledge  regarding  this 
microcontroller  family  played  a  key  role  in  its  selection.  Alt  hough  the  MCU  handles  the  control 
algorithm,  a  lack  of  digital-to-analog  conver  ter  (DAC)  channels  on  commercially  available 
microcontrollers  called  for  implementation  of  an  external  DAC7568  8  -channel  converter.  The 
DAC  chip  was  controlle  d  by  the  MCU  via  a  customized  SPI  interface  in  order  t  o  set  the  DC 
voltage  level  for  each  of  eight  EVP  valves.  Since  the  output  s  of  the  DAC  were  from  0  to  3.3  V, 
supplementary  circuitry  was  required  to  generate  the  10  V  signal  necessary  to  fully  open  the 
EVP  valves.  For  this  purpose,  voltage  gain  wa  s  achieved  through  LF353  operational  amplifiers 
which  were  chosen  for  their  low  cost,  low  drift  a  nd  high  reliability  and  2SD882  tran  sisters  were 
employed  to  provide  t  he  200  mA  current  nece  ssary  to  drive  the  EVP  valves.  With  th  e 
combination  of  the  MCU,  DAC  an  d  amplification  circuitry,  full  control  of  the  EVP  valves  was 
possible. 

Switching  and  Rotating  Circuit:  Currently,  the  overall  system  utilized  a  Spect  rum/Network 
Analyzer,  AC  and  DC  supplies  and  high-power  audio  amplifiers  in  conjunction  with  an  additional 
circuit  which  switched  between  three  detection  coils.  The  switching  circuit  was  com  prised  of  an 
MSP430f2618  which  controls  an  8-channel  multip  lexer.  Simple  digital  swit  ching  between 
Spectrum  Analyzer  inputs  was  made  possible  by  interfacing  the  microcontroller  and  multiplexer, 
with  room  for  additional  sensing  coils  in  the  f  uture.  While  multiplexing  the  input  signals  wa  s 
proven  useful,  the  limitation  of  sam  pling  only  one  coil  at  a  time  was  a  significant  d  rawback  on 
detection  speed.  Additionally,  the  use  of  external  AC  an  d  DC  power  supplies  significantly 
increased  code  complexity  and  system  size. 

Rotation  of  the  detectio  n  platform  was  necessary  to  minimize  the  number  of  det  ection 
coils.  Because  the  response  of  all  strips  from  both  the  front  and  back  were  desired,  a  Portescap 
55M048D2U  stepper  motor  was  implemented  to  rotate  the  platform  between  tests.  The  stepper 
motor  was  controlled  by  the  MCU  through  a  circuit  generated  by  SparkFun  electronics 
(EasyDriver  Stepper  Motor  Driver).  To  avoid  magnetic  interference,  the  power  to  the  motor  was 
disabled  digitally  by  the  EasyDriver  before  measuring  the  sensor  response. 

Microcontroller  Program:  To  maintain  versatility  in  design,  the  MSP430  code  was  I  imited  to  a 
proportional  plus  integ  ral  controller  that  c  ould  receive  commands  from  a  ma  ster  central 
processing  unit  (CPU),  via  standard  RS-232  communication  implemented  on  the  MSP430’  s 
UART  peripheral,  as  well  as  return  load  results.  The  microcontroller  slave  functions  included 
receiving  new  load  data,  using  the  newest  received  load  val  ues,  calibrating  the  load  cell  offsets 
and  sending  the  current  load  state  back  to  the  CPU. 

The  proportional  plus  integral  (PI)  control  fun  ction  was  carried  out  by  modif  ying  the 
DAC7568  output  voltages  via  an  SPI  interface  and  dete  cting  obtained  load  feedback  b  y 
measuring  the  load  cells  connected  to  onboard  ADC  pins.  A  custom  SPI  protocol  was  written  to 
provide  necessary  setup  and  transition  times  as  the  MSP430  peripheral  interface  was  incapable 
of  producing  long  enough  transitions  for  DAC7568  interfacing.  The  control  mechanism  provided 
force  by  varying  the  voltage  level  of  the  EV  P  valve  that  allowed  a  ir  to  exit  th  e  cylinder. 
Increasing  the  voltage  released  more  air  from  the  cylinder  resulting  in  a  decreased  load.  By 
setting  an  initial  voltage  level  at  maximum  outflow  (no  applied  force)  and  adding  the  difference 
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between  the  current  and  the  desired  loads  multiplied  by  so  me  factor  (the  integration  constant), 
an  integral  controller  was  developed. 

Due  to  limitations  of  the  maximum  force  output,  the  inflow  air  level,  a  Iso  controlled  by 
EVP  valves,  was  changed  in  proportion  to  the  desired  load  magnitude.  To  accomplish  this  task, 
calibration  curves  were  generated  which  yielded  the  fastest  settling  tim  es  and  minimal  current 
demands  at  a  desired  load  level.  By  combining  proportional  control  of  the  inflowing  air  rate  with 
integral  control  of  the  outflow  rate,  a  PI  controller  was  generated  which  had  the  capacity  to  apply 
and  measure  loads  ranging  from  0  to  85  Ibf  on  each  of  four  loading  regions. 

In  addition  to  calibrating  the  inflow  air  rate  to  the  desired  load,  consideration  of  the  total 
inflow  rate  for  all  loading  regions  was  necessary.  Because  the  air  supp  ly  was  regulated  for  90 
psi  rather  th  an  flow  rate,  the  flow  rate  varied  according  to  t  he  total  need  of  all  cylinders  at  a 
given  instant.  To  compensate  for  this  inconvenience,  the  inflow  air  rate  was  further  regulated  in 
software  by  increasing  inflow  rate  with  increasing  total  for  ce  application.  This  co  nsideration 
allowed  more  stable  settling  times  across  varied  loading  parameters. 

Computer  Interface:  Although  the  microcontroller  was  designed  to  achieve  optimum  settling 
times  and  stability  across  a  range  of  loading  profiles,  a  computer  interface  was  requ  ired  to  set 
and  update  each  load  as  well  as  recalibrate  and  receive  the  current  load  cell  data.  To  perform 
this  task,  a  Visual  Basic  6.0  (VB6)  interface  was  created  to  communicate  to  the  control  circuitry. 
The  VB6  code  allowed  for  entire  loading  profiles  to  be  actuated  and  recorded  to  a  file  by  simply 
entering  the  settling  time,  starting  load  on  each  cell,  ending  load  on  each  cell  and  total  number 
of  steps  to  collect.  Additionally,  the  code  controlled  the  Spectrum  Analyzer,  AC  and  DC  power 
supplies,  in  order  to  perform  a  DC  biasing  sweep  of  the  loaded  sensor  response,  and  stored  the 
data  to  a  file. 
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Key  Research  Accomplishments 


•  Characterization  3  strip  sensor  grid 

•  Deveiopment  of  superposition  based  aigorithm 

•  Preiiminary  testing  of  etched  sensors 

•  Preiiminary  testing  of  trapezoidai  and  embedded  sensors 

•  Construction  and  testing  of  new  fuiiy  automated  pneumatic  ioader  systems 

•  Fuiiy  automated  simuitaneous  data  coiiection  from  muitipie  coiis 
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Reported  Outcomes 

•  Paper  in  submission  (Appendix  A) 

•  Based  on  the  knowledge  acquired  through  this  work,  we  applied  for,  but  did  not  receive, 
student  fellowships  from: 

o  Michigan  Space  Grant  Consortium 

•  Based  on  the  knowledge  acquired  through  this  work,  we  applied  for,  and  received, 
student  fellowships  from: 

o  National  Defense  Science  and  Engineering  Grant 


14 


Conclusion 


To  date,  continuing  research  has  d  emonstrated  the  efficacy  of  the  sensor  platform  for 
use  on  hard  surfaces  with  preliminary  testing  in  embedded  structures.  At  this  phase,  the  project 
is  heading  towards  testing  which  must  be  pr  eceded  by  experiments  with  non-unif  orm  Metglas 
strips  and  followed  by  more  thorough  testing  with  embe  dded  sensors.  An  automated  fou  r 
cylinder  mechanical  loader  apparatus  was  recently  finishe  d,  with  a  single  force  applicator 
version  already  in  use  for  characte  rizing  a  three  strip  sensor,  and  automated  data  collect  ion 
from  multiple  coils  has  also  been  implemented  in  testing.  Additionally,  preliminary  testing  o  f 
trapezoidal  and  embedded  sensors  was  also  accomplished.  Moreo  ver,  the  e  xcitation  and 
detection  circuits  are  now  ready  for  Integra  tion  and  fu  rther  developments  in  the  use  o  f 
electroplating  and/or  other  processes  to  replace  shearing/etching  sensing  elements  continues. 

These  advances  are  important  due  to  the  lack  of  wireless  long  term  monitoring  systems 
available  to  medical  staff,  not  just  fo  r  lower  limb  prosthetics  but  also  fo  r  other  applications  a  s 
well.  While  the  application  for  this  sensor  syste  m  is  a  lower  limb  prosthesis,  the  technology  can 
easily  be  adapted  to  monitor  other  devices,  such  as  total  knee  arthroplasty  implants.  Moreover, 
the  system  is  expected  to  be  low  cost  and  easy  to  implement  and  produce.  Overall,  this  device 
could  drastically  improve  the  qual  ity  of  care  fo  r  lower  limb  amputees  and  assist  in  the  furthe  r 
development  of  better  prosthesis  in  addition  to  contributing  to  the  understanding  of  the  dynamic 
biomechanical  loads  associated  with  the  stump-socket  interface. 


15 


References 


None 


16 


Appendix  A:  Submitted  Paper  for  Review 


A  Wireless,  Magnetoelastic-based  Sensor  Array  for  Force  Monitoring  on  a  Hard  Surface 

Brandon  D.  Percies,  Andrew  J.  DeRouin,  Thomas  A.  Dienhart,  Ee  Lim  Tan,  Keat  Ghee  Ong* 
Department  of  Biom  edical  Engineering,  Mich  igan  Technological  Einiversity,  Houghton,  MI 
49931,  USA. 

*  Corresponding  author;  kgong@mtu.edu 

Abstract  -  A  force  m  onitoring  system  consisting  of  stress-sensitive  magnetoelastic  strips  for 
remotely  measuring  the  force  profile  across  a  hard  surface  is  described.  Under  the  excitation  of  a 
magnetic  AC  field,  the  m  agnetoelastic  strips  generated  higher-order  harmonic  fields  (magnetic 
AC  fields  at  multiple  frequencies  of  the  excitation  field),  allowing  remote  measurement  of  their 
responses  without  interference  from  the  excitation  field.  Due  to  their  magnetoelastic  properties, 
these  higher-order  harmonic  fields  were  also  dependent  on  the  applied  force  and,  as  a  result, 
variations  in  force/stress  could  be  tracked  vi  a  changes  in  th  e  field  amplitudes.  These  changes 
were  monitored  using  a  detection  system  featuring  a  set  of  m  agnetic  detection  coils,  which 
captured  the  response  of  the  m  agnetoelastic  strips.  To  dem  onstrate  the  functionality  of  this 
sensor  system,  a  three-strip  m  agnetoelastic  sensor  array  was  fa  bricated  on  a  fiat  polycarbonate 
substrate.  The  substrate,  placed  within  a  customized  mechanical  loader,  was  exposed  to  a  variety 
of  force  loading  conditions.  Experim  ental  results  demonstrated  a  proportional  relationship 
between  the  amplitude  of  the  2"*^  order  harmonic  field  and  the  applied  f  orce.  An  algorithm  was 
developed  to  identify  the  m  agnitude  of  the  applied  force.  The  novelty  of  this  system  lies  in  its 
wireless  and  passive  nature,  which  is  ideal  for  applications  in  which  wires  and  internal  power 
sources  are  prohibited  or  discouraged.  Moreover,  the  sensing  component  of  this  system  is  an 
array  of  th  in  magnetoelastic  strips,  allowing  f  or  minimal  modifications  to  existing  structures 
during  implementation. 

Keywords:  Magnetic  higher-order  harm  onic  fields,  m  agnetoelastic,  wireless,  passive,  force 
sensors,  mechanical  loading. 

1.  Introduction 

Eorce  and  stress  are  generally  m  easured  via  strain,  which  is  defined  as  the  change  in 
dimensions  of  an  object  due  to  an  applied  force.  In  practice,  strain  is  often  measured  with  strain 
gauges,  most  of  which  can  be  classified  as  resis  five,  capacitive,  or  vibrational.  Resistive  strain 
gauges  monitor  strain  as  a  function  of  the  change  in  resistance  across  a  conductive  or 
semiconductive  material  when  elastically  def  ormed  [1].  Se  miconductive  materials  exhibit 
piezoresistive  behavior,  resulting  in  a  larger  change  in  the  elect  deal  response  with  applied  force 
compared  to  conductive  materials  [2].  As  a  result,  semiconductive  materials  exhibit  higher  gauge 
factors,  typically  between  50  and  200,  while  c  onductive  materials  experience  less  than  5  [1]. 
These  piezoresistive  strain  gauges  find  common  use  in  different  fields  for  a  variety  of 
applications,  such  as  m  easuring  stress  on  a  knee  prosthesis  [3]  or  stress  m  onitoring  during  an 
electronic  packaging  process  [4].  While  sem  iconductor  strain  gauges  offer  accurate  sensing  on 
small  scales,  their  response  can  vary  with  temperature,  and  they  can  be  difficult  to  manufacture. 

Capacitive  strain  gauges  m  easure  stress/strain  as  a  function  of  the  change  in  the 
capacitance  of  a  sensor.  The  prim  ary  advantages  of  capacitive  strain  g  auges  come  from  their 
capacity  to  operate  in  high  temperature  environments,  minimal  hysteresis,  and  long-term  stability 
[4].  However,  they  have  a  lower  s  ensitivity  compared  to  their  p iezoresistive  counterparts.  Eor 
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example,  thiek  eermet  and  polymer  based  eapaeitive  strain  gauges  demonstrated  similar  linear 
responses  and  hysteresis  to  piezoresistive  sensors  but  with  gauge  faetors  of  only  6  for  the  eermet 
eapaeitor  and  3.5  for  the  polymer  eapaeitor  [5]. 

Vibrating  wire  strain  gauges  function  by  measuring  the  vibrating  frequency  of  a  wire  held 
in  tension  between  two  anchoring  points  and  excited  into  vibration  by  a  magnetic  coil.  Due  to  the 
effects  of  a  pplied  stress/strain,  the  anchorin  g  position  changes,  resulting  in  a  m  easurable 
alteration  in  vibrational  frequency.  [1].  Due  to  their  high  sensitivity,  vibrating  strain  gauges  are 
commonly  used  in  low  strain  structures  such  as  concrete  beams  [6].  As  a  whole,  strain  gauges  in 
this  category  have  the  advantage  o  f  stress  de  tection  in  the  range  of  parts  pe  r  million  and  a  re 
robust,  easily  attachable,  and  ac  curate;  however,  changes  in  tern  perature  can  drastically  affect 
their  performance  [7]. 

Fiber  optic  technology  has  also  been  used  fo  r  stress  monitoring.  The  advantages  of  fiber 
optics  include  electrically  passive  operation,  EMI  immunity,  high  sens itivity,  and  multiplex ing 
capabilities.  Two  m  ain  types  of  fiber  optic  strain  sensor  s  exist:  interferometric  and 
intensiometric.  Interferometric  methods  observe  changes  in  light  passing  through  the  fiber  optic 
cable,  while  intensiometric  sensing  techniques  monitor  changes  in  the  radiant  power  transmitted 
in  a  cable.  A  common  exam  pie  of  an  interfero  metric  sensor  is  the  Fabry-Perot  interferom  eter 
sensor,  which  measures  the  change  in  light  intensity  between  two  mirrors  placed  in  parallel  with 
the  fiber  optic  cab  le.  An  applied  stress  alters  the  distance  between  the  m  irrors  resulting  in  a 
phase  change  in  the  lig  ht.  While  highly  accurate,  this  m  ethod  exhibits  an  inability  to  handle 
periodic  interruptions  by  the  power  supply  [8,9].  Additionally,  fiber  optic  senso  rs  are  prone  to 
damage,  are  affected  by  temperature  elevations,  and  require  that  any  coating  used  to  protect  fiber 
optic  wires  allow  for  proper  transfer  of  force  to  the  sensors  [10]. 

Wireless  capacitive  stress/strain  sensors  were  also  developed  for  monitoring  pressure  as  a 
function  of  change  in  m  ensured  capacitance.  One  such  device  coupled  an  RE  transceiver  to  a 
custom  capacitor  for  pressure  monitoring  [11].  The  device  was  comprised  of  a  pressure  sensitive 
cavity,  fabricated  through  silico  n  fusion  bonding  of  two  silicon  wafe  rs  used  to  seal  the  cavity, 
and  a  flexib  le  and  a  st  ationary  electrode,  form  ing  a  capacito  r.  Pressure  applied  to  the  device 
deflected  the  flexible  electrode  toward  the  fixe  d  electrode,  thus  altering  the  capacitance  of  the 
sensor.  Using  a  RE  transceiver,  inform  ation  was  wirelessly  collected  and  sent  for  analysis  [11]. 
Another  wireless  capacitive  sen  sor,  known  as  th  e  SmartPill,  incorporated  pressure,  pH,  and 
temperature  sensors  for  m  onitoring  gastrointestinal  tract  conditions  [12].  Unfortunately,  these 
systems  are  limited  in  their  size  by  the  necessity  for  onboard  electronics  and,  in  the  case  of  an 
active  sensor,  may  eventually  require  the  device  to  be  removed  in  order  to  replace  the  power 
supply. 

A  simpler  version  of  wireless  capacitive  stress/strain  sens  ors  was  realized  by 
incorporating  an  inductive-capaci  five  (EC)  tank  circuit  into  a  sensor.  In  a  wireless,  passive 
stress/strain  EC  sensor,  the  capac  itive  stress/strain  element  connected  to  an  indu  ctor,  which 
remotely  conveyed  stress/strain  inform  ation  as  a  change  in  the  resonant  frequency  of  the  tank 
circuit.  Among  the  applications  of  this  type  of  sensor  are  monitoring  st  ent  integrity  after  an 

endovascular  repair  procedure  [13]  and  measuring  pressure  in  automobile  tires  [14]. 

Another  class  of  passive  st  ress/strain  sensors  is  base  d  on  am  orphous  magnetoelastic 
materials.  When  exposed  to  a  time  varying  AC  field,  magnetoelastic  materials  vibrate  due  to  the 
magnetoelastic  effect.  The  m  agnetoelastic  effect  also  causes  the  vibrating  m  agnetoelastic 
material  to  generate  a  magnetic  flux  that  reaches  a  peak  at  its  m  echanical  resonant  frequency 
[15,16].  When  an  internal  stress  is  applied,  the  resonant  frequenc  y  of  the  strip  becom  es  stress 
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dependent,  and,  as  an  exam  pie,  has  been  used  to  determine  atmospheric  pressure  [15,16].  This 
sensor  is  not  only  passive  and  wireless,  but  also  low  cost  and  long  las  ting.  However,  the  strip 
itself  must  be  stressed  (by  bending)  to  act  as  a  se  nsor  and  there  is  no  relia  ble  way  to  control  its 
sensitivity. 

A  strain  sensor  was  developed  using  magnetoelastic  materials  by  Kouzoudis  and 
Mouzakis  [17].  A  Metglas  2826MB  ribbon  was  attached  on  an  epoxy  resin  slab  and  exposed  to 
vibrations  of  varying  amplitudes  and  frequenc  ies.  Due  to  its  m  agnetoelastic  property,  the 
vibration  of  the  magnetoelastic  material  caused  a  change  in  its  magnetization  stages,  which  was 
remotely  picked  up  by  a  nearby  de  tection  coil.  The  sensor  was  demonstrated  to  have  a  strain 
gauge  factor  of  1 1,700  at  a  vibration  frequency  of  150  Hz. 

Due  to  magnetic  softness,  amorphous  magnetoelastic  materials  also  generate  higher  order 
harmonic  fields  (magnetic  fields  at  multiple  fre  quencies  of  the  excitatio  n  field)  when  under  the 
excitation  of  a  low  frequency  AC  magnetic  field  [18].  To  visu  alize  the  higher-order  harm  onic 
fields,  the  magnetoelastic  material  is  generally  excited  by  a  steady  AC  magnetic  field  along  with 
a  sweeping  DC  biasing  field.  The  biasing  field  alters  the  magnitude  of  the  higher  order  harmonic 
fields  and  thus  produces  a  di  stinct  pattern  as  shown  in  Figure  1.  In  the  previous  work,  it  was 
shown  that  the  amplitude  of  the  higher-order  harmonic  fields  increased  when  force  was  applied 
along  the  length  of  a  magnetoelastic  strip  [19]. 

This  paper  describes  a  force  m  onitoring  system  that  tracked  the  changes  in  the  2"‘*  order 
harmonic  amplitudes  of  an  array  of  m  agnetoelastic  strips.  By  using  an  array  of  magnetoelastic 
strips,  the  new  system  is  able  to  m  onitor  not  only  the  total  contact  force  on  the  surface,  but  can 
also  determine  the  force  distribution  on  the  surface.  In  contrast  to  the  previous  work  [18],  fore  e 
was  applied  directly  onto  the  surface  instead  of  along  the  lengt  h  of  the  m  agnetoelastic  strip. 
Similar  to  the  previous  work,  however,  the  ap  plication  of  force  on  th  e  magnetoelastic  strip 
surface  resulted  in  an  increase  in  the  2“‘*  order  harmonic  field  (see  Figure  1). 

By  capturing  the  respo  nse  of  each  strip,  the  system  could  identify  the  position  and 
magnitude  of  the  applied  force.  This  sensing  sy  stem  was  not  only  wireless  and  passive,  but  also 
simple  to  implement  since  it  was  nothing  m  ore  than  an  array  of  m  agnetoelastic  strips  directly 
applied  onto  the  desired  surface.  The  process  re  duces  cost  and  will  allow  for  the  production  of 
long  lasting  sensors  for  a  variety  of  applications. 

It  is  wor  th  noting  that  the  presented  sensor  differs  from  other  devices  based  on 
magnetoelastic  materials,  such  as  those  described  by  Kouzoudis  and  Mouzakis  [17],  in  terms  of 
operating  principle,  ideal  application,  and  strengths  and  weaknesses.  For  instance,  to  obtain  good 
sensitivity,  Kouzoudis’  sensor  system  required  the  substrate  to  vibr  ate  and  the  sensor  sensitivity 
was  also  proportional  to  the  vibrational  frequency.  The  described  sensor,  on  the  other  hand,  does 
not  require  the  substrate  to  vibrate,  and  its  sensitivity  is  largely  related  to  the  magnetoelasticity  of 
the  material.  Moreover,  unlike  Kouzoudis’  system  ,  the  described  sensor  uses  the  higher-order 
harmonic  signals  from  the  m  aterial  to  track  pressure  and  st  ress.  The  use  of  higher-order 
harmonic  signals  can  significantly  rem  ove  the  background  excitation  signal,  thus  increasing  the 
signal  to  noise  ratio.  Furthermore,  the  major  applications  for  both  technologies  are  different  due 
to  the  differences  in  their  operating  principle.  Thus,  the  described  sensor  is  better  suited  for  stress 
and  pressure  monitoring  while  Kouzoudis’  sensor  system  is  ideal  for  strain  monitoring. 

2,  Experiments 
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2.1,  Sensor  Fabrication  and  Experimental  Setup 

The  sensor,  illustrated  in  Figure  2,  was  eut  from  a  block  of  polycarbonate  material  and 
measured  47  mm  x  56.5  mm  x  12.2  mm.  Metglas  2826MB  ribbon  (Fe  4oNi38Mo4Bi8),  purchased 
from  Metglas  Inc,  Conway,  SC,  USA,  was  used  as  the  magnetoelastic  stress  sensing  material  due 
to  its  large  m  agnetostriction  (>12  ppm  ),  high  permeability  (>50,000),  and  low  m  agnetic 
coercivity.  Three  sensing  strips  were  sheared  from  a  26  pm  thick  Metglas  ribbon  to  50  mm  x  5 
mm,  and  were  adhered  to  the  upper  surface  of  t  he  block  using  cellulose  adhesive  tape  (50  pm 
thick).  Adhesive  tape  was  found  to  introduce  less  internal  stre  ss  than  glue  or  epoxy,  thus 
preserving  sensor-to-sensor  reproducibility.  For  convenience,  these  sensing  strips  were  labeled 
Strip  1,  2,  and  3  respectively. 

An  automated  pneumatic  mechanical  loader  (s  ee  Figure  3)  was  constructed  to  apply 
controllable  force  to  each  strip  of  the  sensor.  Figure  3  also  illustrates  the  control  elements  of  the 
automated  system.  A  manual  air  control  valve  allowed  for  adjustments  to  the  overall  air  flow  into 
the  apparatus.  Following  the  m  ain  valve,  the  air  was  split  between  an  E  VP  series  Proportional 
Control  Valve  (Clippard  Instrum  ent  Laboratory,  Inc.  7390  Colerain  Ave.,  Cincinnati,  Ohio 
45239)  and  a  pneumatic  piston.  Opening  or  closi  ng  the  EVP  valve  by  incr  easing  or  decreasing 
the  voltage  from  a  connected  Kepco  lOV  Programmable  Power  Supply  altered  the  volume  of  air 
flowing  to  the  pneum  atic  piston.  This  in  turn  affected  the  force  applied  by  the  piston.  Applied 
force  was  monitored  u  sing  a  Measurem  ent  Specialties  PC23  Compression  Load  Cell,  placed 
between  the  load  applicator  plate  and  the  pi  ston.  An  Extech  382202  DC  Power  Supply  provided 
the  input  voltage  to  the  load  cel  1  and  applied  force  data  was  m  onitored  and  transferred  to  a  PC 
using  a  Hewlett  Packard  3478 A  Multim  eter.  A  custom  Visual  Basic  program  calculated  the 
difference  between  the  actual  applied  force  and  the  desired  force  which  was  then  altered  using  an 
integral  control  algorithm  to  produce  an  appropri  ate  change  in  voltage  to  be  sent  to  the  EVP 
valve.  The  system  provided  repeatable  loading  with  a  range  of  0-355.86  N  and  was  capable  of 
adjusting  the  force  by  4.44  N  ±  1.1 1  N.  In  addition  to  autom  ated  force  loading,  the  system  and 
Visual  Basic  program  allowed  for  simultaneous  collection  of  sensor  responses  from  the  detection 
coils.  It  is  worth  noting  that  the  range  and  wei  ght  increment  of  the  system  were  set  for  the 
purposes  of  this  testing  by  adjust  ing  the  air  pressure  into  the  sy  stem,  the  air  flow  at  the  m  ain 
regulator,  the  coefficient  used  in  calculating  the  change  in  voltage  to  the  EVP  valve,  and  the  time 
allowed  between  weight  increments  for  the  system  to  settle. 

Rectangular  detection  coils  consisted  of  a  functioning  coil  connected  to  an  oppositely 
wound  compensating  coil  in  series,  both  m  ade  of  100  turns  of  36  gage  copper  wire.  The 
individual  coils  measured  12.0  mm  x  17.6  mm  x  4.3  mm.  Eor  convenience,  the  detection  coils 
were  labeled  Coil  1,  2,  and  3  acco  rding  to  the  strip  bein  g  measured  (see  Figure  2),  and  the 
compensating  coils  we  re  similarly  labeled  Compensating  Coil  1,  2,  and  3.  During  the 
experiments.  Strip  1,  2,  and  3  were  aligned  to  the  centers  of  Coil  1,  2,  and  3,  respectively. 

The  sensor  was  secured  within  the  mechanical  loader  and  the  a  pparatus  was  positioned 
directly  in  front  of  the  detection  and  excitation  coils.  As  illustr  ated  in  Figure  4,  the  detection 
coils  were  connected  to  an  Agilent  spectrum  /network  analyzer  4396B  to  capture  the  signal  for 
the  PC  (through  a  custom  Visual  Basic  program  and  GPIB  interface)  for  further  analysis.  The 
excitation  coils  consisted  of  two  superim  posed  50  turn  18-gauge  coils  (28  cm  in  diameter)  that 
provided  the  AC  and  DC  excitati  on  fields.  One  coil  was  con  nected  to  an  AC  function  generator 
(Fluke  271  10  MHz)  and  an  a  mplifier  (Tapco  J 1400),  while  the  other  coil  (DC)  connected  to  a 
Kepco  MBT  36-10  MT  power  supply.  In  all  tests  the  AC  field  was  150  A/m,  200  Hz  and  the  DC 
field  was  0-250  A/m. 


20 


2,2,  Experimental  Procedure 

The  response  of  eaeh  m  agnetoelastic  strip  wa  s  colleeted  while  the  total  app  lied  load 
increased  from  0  to  266.89  N  and  then  decreased  to  0  N  at  weight  increments  of  22.24  N.  During 
the  experiment,  data  was  collected  sim  ultaneously  from  all  coils  following  a  short  period  to 
allow  the  piston  to  reach  the  desired  load  and  stabilize  there.  The  collected  data  was  then  zeroed 
to  a  common  starting  point  by  subtracting  a  zero  load  value  obtained  from  the  first  data  point. 
This  procedure  was  repeated  for  a  variety  of  lo  ading  conditions  created  by  placing  rubber  inserts 
over  specified  strips,  thus  produc  ing  conditions  where  in  som  e  strips  were  loaded  and  others 
were  not. 

In  addition  to  changing  with  the  application  of  force,  the  response  of  the  sensing  strips 
also  varied  as  a  function  of  relative  location  from  the  detection  coils.  To  investigate  the  effect  of 
changing  sensor  location,  strip  responses  were  measured  while  m  oving  the  m  echanical  loader 
incrementally  on  the  x,  y,  and  z  axes  of  a  rectangular  coordina  te  system.  The  origin  of  the 
coordinate  system  was  defined  as  the  exact  center  between  Coil  2  and  Compensating  Coil  2  (see 
Figure  3  and  Figure  4). 


2.3.  Stress  Sensing  Theory 

The  pressure  sensitivity  of  the  magnetically  soft  magnetoelastic  material  can  be  explained 
by  the  m  agnetic  susceptibility  ( ;f),  which  is  the  ra  tio  of  saturation  magnetization  to 
anisotropy  field  (T4)  expressed  as  [20]; 

x  =  mjh,  (1) 

The  anisotropy  field  of  a  m  agnetic  material  can  be  related  to  the  te  nsile  stress  along  the 
magnetization  direction  as  [16]; 

H,=H,,-U,c7^IM^  (2) 

where  Hko  is  the  anisotropy  field  at  zero  stress,  Xs  is  the  saturation  magnetostriction  of  the 
material,  and  oy  is  the  tensile  stress  along  the  m  agnetization  direction,  which  is  also  along  the 
length  of  the  sensor. 

Eq.  (2)  describes  the  change  in  anisotropy  field  due  to  the  tensile  stress  along  the  sensor’s 
length;  however,  for  this  particular  application,  force  was  loaded  on  the  dominant  surface  of  the 
ribbon  shape  sensor  (along  the  z-direction  shown  in  Figure  2).  Therefore,  the  transverse  stress  on 
the  sensor  surface  (  z-direction)  was  related  to  the  tens  ile  stress  along  the  sensor  length  (  y- 
direction)  using  the  Poisson’s  ratio  o  as; 

ct^=2(tJv  (3) 


Note  that  a  scaling  factor  of  two  was  added  in  Eq.  (3)  to  com  pensate  for  the  fact  that  only  one 
side  of  the  sensor  was  being  stressed. 

As  shown  in  Eq.  (2),  increasing  stress  d  ecreases  the  anisotropy  field  of  a  m  agnetic 
material,  assuming  the  anisotropy  energy  and  the  saturation  magnetization  stay  constant.  The 
change  in  anisotropic  field  has  a  direct  impact  on  the  measured  signal  amplitude  of  the  n-th  order 
magnetic  harmonic  field  {An),  in  Volts,  which  can  be  described  by  the  equation  [17]; 
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where  E  is  a  variable  that  accounts  for  sensor-coil  coupling,  Bs  is  the  satura  tion  induction  flux 
(for  ferromagnetic  materials,  B^  »  ),  o)  is  the  radian  frequency  of  the  fundamental  order,  hac  is 


21 


the  AC  excitation  field,  and  Hdc  is  the  DC  biasing  field.  Under  an  externally  applied  load,  the 
sensor  material  deforms  and  generates  an  internal  stress.  Eqs.  (1)  -  (4)  indicate  that  stress  causes 
a  change  in  m  agnetic  anisotropy,  which  alters  m  agnetization  and  changes  the  higher-order 
magnetic  fields  allowing  for  remote  detection  of  pressure. 

As  indicated  in  Eqs.  (1)  and  (2),  x  is  inversely  proportional  to  E4,  which  decreases 
linearly  with  increasing  cr^.  As  a  result,  the  susceptibility  of  the  material  is  expected  to  show  an 
exponential  pattern  with  increasing  stress  and  eventually  converge  on  an  asymptote  at  infinity  as 
Hk  approaches  zero.  While  theoretically  this  represents  an  infinite  increase  in  susceptibility,  and 
thus  an  infinite  increase  in  the  2  order  harmonic  amplitude,  realistically  the  change  in  the 
susceptibility  and  the  2"‘*  order  harmonic  amplitude  will  more  likely  resemble  an  upper  bounded 
decaying  exponential  curve  since  the  susceptibility  will  experience  a  much  slower  change  when 
the  anisotropy  has  reached  a  near  zero  state. 

2.4.  Determination  of  the  Force  Loading 

An  algorithm  was  developed  to  determine  the  force  loading  on  each  magnetoelastic  strip 
based  on  the  m  easured  2”‘*  order  harmonic  amplitudes.  Since  the  measured  2“  order  harmonic 
amplitudes  of  the  m  agnetoelastic  strips  we  re  expected  to  f  ollow  an  upper  bounded  decaying 
exponential  function  with  increasing  stress,  upper  bounded  decaying  ex  ponential  curves  were 
used  to  fit  the  m  easured  data.  Due  to  the  clos  e  proximity  of  the  m  agnetoelastic  strips  to  one 
another,  the  stress  resp  onse  of  each  m  agnetoelastic  strip  was  also  dependent  on  the  loading 
conditions  of  its  neighboring  strips.  This  resu  It  is  expected,  based  upon  the  sim  pie  fact  that 
neighboring  magnetic  fields  will  interact  with  one  another.  Therefore,  assuming  the  strip-to-strip 
cross  interference  is  cumulative  at  a  g  iven  detection  coil,  the  measured  2"^*  order  harm  onic 
amplitude  of  magnetoelastic  Strip  i  (measured  by  the  detection  coil)  was  represented  by  the 
summation  of  the  responses  of  all  strips  as: 

(5) 

./=1 

where  St  is  the  m  easured  2"‘^  order  harm  onic  amplitude  at  Coil  i.  Ay  is  the  pea  k  amplitude 
measured  at  Coil  i  when  only  strip  j  is  at  the  maximum  loading,  ay  is  the  decay  co  efficient  (at 
Coil  i  when  strip  j  is  loaded),  and  f  is  the  applied  force  at  Strip  j. 

Eq.  (5)  consists  of  three  upper  bounded  de  caying  exponential  equations  with  three 
unknowns.  To  solve  Eq.  (5),  a  sim  pie  iterative  method  was  developed  to  identify  f  for  a  set  of 
given  Si.  Starting  with  a  set  of  es  timated  force  load  ing,  the  iterativ  e  process  determines  the 
difference  {di)  between  the  calculated  Si  from  Eq.  (5)  and  the  measured  signal  Si  as: 

3 

5^  =  S- {Calculated) - S. {Measured)  =  ) ~ Sj {Measured)  (6) 

7=1 

A  zero  5  for  all  strips  indicates  the  correct  input  for  the  force  loadings;  in  contrast,  a  non-zero  5 
indicates  there  is  an  error  in  the  calculated  force  load  ings.  The  erron  eous  force  loading  gj  is 
determined  from  di  as: 

S,=A.(l-e’'A  (7) 

The  new  estimated  force  loading  is  determined  by  subtracting  the  current  force  loading  to 
the  erroneous  force  loading  gf. 
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Eq.  (6)  and  (7)  are  itera  tively  solved  until  /4' 
iterative  solution. 


-fj<£,  where  £  is  the  aeeeptab  le  error  of  the 


3.  Results  and  Discussion 

Figure  5  plots  the  signal  recorded  by  Coil  1,  2,  and  3  when  rubber  inserts  were  placed  on 
Strip  1  (Figure  5a),  2  (Figure  5b),  or  3  (  Figure  5c),  respectively.  These  results  indicate  that 
when  the  magnetoelastic  strips  were  under  direct  force  loading,  the  si  gnal  recorded  by  their 
corresponding  coils  in  creased  following  an  upper  bound  ed  decaying  exponential  function  : 
a(\-e  '‘),  confirming  the  theoretical  behavior  of  the  strips  described  previously.  Also  noticeable 
is  that  in  Figure  5a,  the  stress  response  m  ensured  at  Coil  2  was  not  zero,  but  instead  was  about 
25%  of  the  stress  response  measured  at  Coil  1,  even  though  Strip  2  was  not  loaded.  This  is  due  to 
interference  from  Strip  1.  In  contrast.  Coil  3  measured  a  zero  response  because  the  interference 
from  Strip  1  was  shielded  by  Strip  2  before  reaching  Strip  3.  Figure  5c  shows  a  similar  response 
as  Figure  5a  since  the  sensor  was  symme  trical;  however,  the  response  at  Coil  3  was  slightly 
lower  compared  to  Coil  1  due  to  minor  differences  in  physical  dimensions  between  Coil  1  and  3. 

Similarly,  Figure  5b  indicates  that  Strip  2  interfered  w  ith  measurements  at  Coil  1  and  3. 
Ideally,  the  measurements  at  Co  il  1  and  3  should  be  identical  but  the  response  at  Coil  3  was 
slightly  smaller  than  Coil  1  due  to  d  ifferences  in  coil  dimensions.  Also,  compared  to  Figure  5a, 
the  change  in  harmonic  amplitude  was  higher  when  Strip  2  was  loaded  since  Strip  2’s  response 
was  interfered  with  by  both  Strip  1  and  3,  while  in  Figure  5a  only  Strip  1  interfered  significantly 
with  Strip  2. 

After  determining  the  coefficients  Ay  and  Uy,  the  performance  of  the  iteration  process  was 
examined.  To  prevent  m  easurement  errors  from  affecting  the  iteration  process,  all  input 
measurements  Si  were  calculated  using  Eq.  (5).  Figure  6  plots  the  estimated  force  loadings  for 
all  strips  at  each  iteratio  n  step  (initial  force  loadings  were  set  to  zero).  It  was  fou  nd  that  the 
performance  of  the  iteration  process  degraded  with  increasing  force  loading  and  that  in  the  worst 
case  scenario  (156.58  N  loading),  the  error  (  s)  was  0.35%  after  1000  iterations  and  0.0075% 
after  2000  iterations. 

To  determine  the  performance  of  the  whole  sensor  system,  all  magnetoelastic  strips  were 
loaded  and  the  responses  from  Coil  1,  2,  and  3  were  m  ensured  simultaneously.  The 
measurements  were  then  fed  into  Eqs.  (6)  and  (8)  to  ite  ratively  solve  for  the  force  loading  on 
each  strip.  Figure  7  plots  the  absolute  percentage  error  between  the  actual  and  calcu  lated  forces 
on  all  strips  at  different  force  loading  condition  s.  Although  the  iteration  process  has  an  error  of 
only  0.0075%  after  2000  it  erations,  due  to  uncertainties  in  the  m  ensured  data  and  other 
experimental  errors,  the  calculated  force  has  a  1 0%  error  when  compared  to  the  measured  force. 
One  source  of  error  was  from  the  current  method  of  force  application.  Rubber  inserts  were  used 
to  distribute  load  on  the  strips.  When  the  rubber  inserts  were  placed  on  different  strips,  the  force 
distribution  was  assumed  to  be  equally  distributed  among  the  strips.  In  practice,  some  strips  may 
experience  more  force  than  the  oth  ers  due  to  slight  imbalances  of  the  loading  plate.  As  a  result, 
the  measured  force  on  each  strip  m  ight  vary  slightly  from  the  actual  load  on  the  strip.  Another 
source  of  error  was  the  exponential  na  ture  of  the  data.  As  shown  in  Figure  5,  the  sensitivity  of 
the  sensor  decreased  with  increasing  force  loading  as  the  upper  bounded  decaying  exponential 
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curve  saturated.  As  a  result,  the  aeeuraey  of  the  sensor  deereased,  espeeially  in  high  load  regions, 
leading  to  larger  errors. 

In  addition  to  eharaeterizing  the  sens  or  response  under  loading  and  analyzing  the 
developed  algorithm,  the  effeets  of  movement  on  sensor  response  were  also  analyzed.  The  signal 
at  Coil  2  was  reeorded  and  presented  in  Figure  8  when  the  sensor  was  inerem  entally  moved  on 
the  X,  y,  and  z  axes.  While  moving  along  the  x  axis,  the  measured  amplitude  inereased  when  eaeh 
sensing  strip  was  eloser  to  the  een  ter  of  the  eoil  and  then  deereased  as  that  eenter  p  assed.  This 
result  was  observed  three  tim  es  (eorresponding  to  the  3  strips)  followed  by  a  sharp  deerease 
toward  zero.  A  si  milar  result  was  seen  along  the  z  axis;  however,  in  this  ease  the  rising  and 
falling  of  the  response  o  eeurred  as  Strip  2  m  oved  within  Coil  2  and  C  ompensating  Coil  2.  As 
expeeted,  they  axis  testing  demonstrated  a  gradual  signal  deerease  with  distanee. 

From  Figure  8,  it  is  evident  that  the  loeation  of  the  deteetion  eoils  played  a  eritieal  role  in 
the  aeeuraey  of  the  system.  In  the  experiment,  the  loeation  of  the  sensor  was  fixed  with  respeet  to 
the  deteetion/exeitation  eoils.  However,  in  praetieal  use,  the  user  may  not  be  able  to  position  the 
deteetion  eoils  aecurately.  To  ensure  integrity  of  the  measurement,  the  user  will  be  able  to  move 
the  eoils  around  the  sensor  until  the  system  pieks  up  the  maximum  signal,  indieating  the  eorreet 
orientation  of  the  sensor.  Altem  atively,  it  is  pos  sible  to  plaee  a  ealib  ration  sensing  strip,  whieh 
ean  be  parallel  to  the  stress-responsive  sensing  strips  but  at  a  position  that  is  insulated  from  foree 
loading,  sueh  that  all  m  easurement  data  is  ealibrated  from  the  ealibration  sensing  strip  to 
eliminate  the  loeation  effeet. 

4.  Conclusion 

The  fabrieation  and  testing  of  a  wireless  passive  sensor  system  for  monitoring  the  applied 
foree  on  a  hard  surfaee  has  been  presented.  Th  e  sensor  was  tested  from  0  to  226.  89  N  with 

experimental  data  demonstrating  an  exponential  inerease  in  the  2  order  harmonie  amplitude  of 
magnetoelastie  sensing  strips  as  pressure  inereased  .  It  was  dem  onstrated  that  neighboring  strips 
have  an  interferenee  effeet  wh  ieh  eontributes  to  the  overall  se  nsor  reading  at  a  given  eoil. 
Additionally,  a  simple  iterative  algorithm  was  developed  to  determine  the  applied  f  oree  on  all 
sensing  strips  by  examining  the  signals  eaptured  by  the  deteetion  eoils. 

Future  works  inelude  the  design  and  fabrieation  of  a  more  eomplieated  sensor,  whieh  will 
inelude  strips  forming  a  sensing  grid,  and  a  m  ore  sophistieated  algorithm  to  aeeommodate  the 
more  eomplex  sensor  strueture.  In  addition,  a  new  proeess,  sueh  as  eleetroplating  or  sereen 
printing,  will  be  developed  f  or  fabrieation  of  the  sensin  g  strips  to  reduee  sens  or-to-sensor 
variability. 
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Figure  Captions 

Figure  1 .  The  (a)  2  order  harm  onic  field  and  (b)  3  order  harm  onic  field  measured  as  a 
function  of  an  applied  DC  field  with  and  without  an  applied  force. 

Figure  2.  The  sensor  was  com  prised  of  three  magnetoelastic  sensing  strips  affixed  to  a 
polycarbonate  substrate. 

Figure  3.  The  mechanical  loader  housed  the  sensor  and  allowed  for  incremental  loading  of  the 
sensing  strips.  The  total  applied  force  was  measured  with  a  load  cell. 

Figure  4.  The  full  experimental  setup  illustrating  the  excitation  coils  and  the  detection  co  ils. 
During  the  experim  ents,  the  m  echanical  loader  (not  shown  here)  and  the  test 
substrates  were  placed  directly  in  front  of  the  detection  coils  so  Strip  1,2,  and  3  were 
directly  aligned  to  Coil  1,  2,  and  3,  respectively. 

Figure  5.  Changes  in  the  2"‘*  order  harmonic  amplitude  captured  by  Coil  1,  2,  and  3  when  onl  y 
Strip  1  (a),  2  (b),  or  3  (c)  was  loaded,  respectively. 

Figure  6.  Force  loading  on  the  m  agnetoelastic  strips  estimated  by  the  iterative  process  as  a 
function  of  iteration  step.  The  numbers  in  the  figure  are  the  expected  forces. 

Figure  7.  The  percentage  error  between  the  actual  and  calculated  forces  on  all  strips  at  different 
force  loading  conditions. 

Figure  8.  Response  of  sensor  when  incrementally  moved  along  the  x,  y,  and  z-axes. 
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Figure  lb. 
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Figure  4. 
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Figure  5b. 
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